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ABSTRACT 

The 11.3 /iui emission feature is a prominent member of the family of uniden- 
tihed infrared emission (UIE) bands and is frequently attributed to out-of-plane 
bending modes of polycyclic aromatic hydrocarbon (PAH) molecules. We have 
performed quantum mechanical calculations of 60 neutral PAH molecules and 
found that it is difficult to reconcile the observed astronomical feature with any 
or a mix of these PAH molecules. We have further analyzed the fitting of spectra 
of several astronomical objects by the NASA PAH database program and found 
that reasonable fittings to the observed spectra are only possible by including 
significant contributions from oxygen and/or magnesium containing molecules in 
the mix. A mixed of pure PAH molecules, even including units of different sizes, 
geometry and charged states, is unable to fit the astronomical spectra. Prelimi¬ 
nary theoretical results on the vibrational spectra of simple molecules with mixed 
aromatic/aliphatic structures show that these structures have consistent bundles 
of vibrational modes and could be viable carriers of the UIE bands. 

Subject headings: infrared: ISM - ISM:hnes and bands — ISM: molecules — 

ISM: planetary nebulae: general 
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Introduction 


The unidentified infrared emission (UIE) bands at 3.3, 6.2, 7.7, 8.6, and 11.3, and 
12.7 /um are a family of infrared emission bands which are widely observed in planetary neb¬ 
ulae, reflection nebulae, Hll regions, diffuse interstellar clouds, and galaxies. In star forming 
galaxies, UIE ban ds contribute to up to 20% of the total infrared output of the galaxies 
(jSmith et al.l 1200711 . The carrier of the UIE features therefore must be a major constituent 
of interstellar matter and their spectral characteristics can serve as a pot ential probes for 
the chemical proces ses and physical conditions of interstellar environments (iJoblin fc Tielens 


fr K 

20111: iPeetersI 1201311 . 


The 11.3/Ltm feature is a prominent me mber of the UIE family and was among the 
first UIE bands detected flGillett et aI.Nl973[l . The detection of this feature in the plane¬ 
tary nebulae NGC 7027, BD-|-30°3639, and NGG 6572 came as a surprise as only atomic 
fine-structure lines were expected to be present in these objects. The strength of the 11.3 
/rm feature was found not to be related to the excitation conditions of the nebulae an d 
its origin was first attributed to mineral carbonates such as MgGOs flGillett et al.Nl973l) . 
This identificati on was soon discard ed when the expected 7 .0 pm feature from carbonates 
was not found flRussell et ahl Il977bfl . iRussell et al.l f)l977all found the 11.3/rm feature to 
be closely correlated with the 3.3 /xm feature, suggesting a common origin for the two fea¬ 
tures. The current accepted explanation is that these two featu res arise from the G-H 
vibration modes of aromatic compounds flDulev fc WilliamsI Il98ll) . Many aromatic com¬ 
pounds are known by laboratory spectroscopy to show absorption bands around 3 /im due 
to G-H stretch and around 12 pm due to G-H out-of-plane (OOP) bending modes. The 
fact that the 11.3 pm feature is found almost exclusively in planetary nebulae with high 
G to O ratio (> 2), and sometime s associated with the faint SiG feature, suggests carbon 
grains as the carrier flBarlowl 1198311 . Since then, various carbonaceous materials containing 
an aromatic component have been suggested as the UIE carrier. The li st of candidates in¬ 
cludes hydrogenate d amorphous carb on IHAG. iDulev fc WilliamsI 1198311 . quenched carbon 
composites (QCC, 


flheger fc Pugetlll984 : 
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1987h. poly ey e. 
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19891) . petroleum fractions fiGataldo. &: Kehevan 
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ic aromatic hydroca rbon ( PAH) molecules 
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matic/aliphatic organic nanoparticles fMAONs. iKwok &: Zhang 


2011 


coal fjPapoular et ah 


2013 h and mixed aro- 


20131) . Among these 


proposals, the PAH hypothesis has captured the most attention and is widely believed as 
the carrier of the UIE features by the astronomical community. 

Assuming a PAH origin, the 11.3 pm feature ha s been frequently used as a tool to probe 


the physical co nditions of galaxies and AGNs (e.g., ISmith et al.l 120071: iGalliano et ah 


2008 


Wu et al.ll2010l) . The 11.3 pm feature is used to trace the neutral and/or large PAHs, and 
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the strength ratios between the 11.3/rm featnre and shorter wavelength (snch as 6.2, 7.7, 
and 8.6 /im) UIE bands have been nsed to estimate the properties of the radiation helds. 
The physical validity of these stndies, however, depends on a correct interpretation of the 
origin of the 11.3 pm featnre. 


Recent advances in quantnm-chemical methodology and compntational power provide 
the accuracy needed to quantitativel y study the UIE phenornenon. The NASA Am es PAH 
IR Spectroscopic Database fPAHdb: iBauschlicher et al.ll2010l: iBoersma et al.ll2014l) collects 
theoretical infrared spectra of 700 PAH molecules. Using the density functional theory 
(DFT), we have cal culated the IR spec tra of a series of molecules with mixed aliphatic- 
aromatic structures fISadiadi et al.ll2015l) . In this study, we reexamine the 11.3 pm feature 
utilizing the PAHdb data and new quantum chemistry calculations. 


2. Observations of the 11.3 pm feature in different astronomical sources 


The observed peak position of the 11.3 pm feature in astronomical sources is well de- 
hned and d oes not vary rnuch i n wavelength. The 11.3 pm feature is also detected in 
absorption flBregman et al.N2000l) . and its wavelengths and prohles closely resemble those 
seen in emission. The 11.3 pm feature has a distinctive asymmetric profile, having a steep 
decline in the short wavelengt h side and a gradual extended wing in the long wavelength side 
(Ivan Diedenhoven et al.ll2004j) . Such asymmetric profiles are difficult to explain by gas-phase 
molecular emissions. Explanations offered include having PAH molecules of different mass 
ranges, with the high-mass PAHs r esponsible for the short -wavelength side and low-mass 
PAHs for the long-wavelength side (jCandian Sz Sarrd l2015h . Furthermore, the peak posi¬ 
tions of the C-H OOP bending modes of PAH molecules can be quite different du e to mode 
coupling, molecular structure and charge states of the molecules flHonv et al.ll200ll) . Because 
of the large wavelength variations of the OOP modes of PAH molecules, it is difficult to as¬ 
sign or match the observed astronomical feature to specihc PAH molecules. One suggestion 
of a possible s olution to this proble m is to incorporate internal hydrogen in carbonaceous 
microparticles ( Balm fc Krotdll9901) . Ring deformation vibrational mode of small carbona¬ 
ceous molecules such as ethylene oxide (c-C-;iH 40 ') has also been proposed to explain the 
narrow 11.3 pm feature (IBernstein fc Lvnchl 120091 1. 


Like the other UIE bands, the 11.3 pm feat ure is not seen in A GB stars but only 
emerges during the post-AGB phase of evolution. iGoebel et ahl fjl995l) did not detect this 
feature in 718 carbon stars, except a feature at 11.9 pm in several sources. The situation 
is further complicate d by the blending of a SiG feature (“class 5” spectra according to 
Matsuura et ahl 120141) . A sample study of Herbig Ae/Be stars shows that the central wave- 



































4 


length of the 11.3 am fe ature decreases with increasing effective temperature of the central 
stars flKeller et al.l 120081) . suggesting that the chemical structure of its carrier is being pro¬ 
cessed by stellar radiation. In the spectral classihcation scheme of iPeeters et al.l fl2002l) . the 
11.3 /im feature occurs at longer wavelength in Class C sources in comparison to those in 
class A and B objects. Since Class C sources usually show aliphatic feature (e.g., at 3.4 /rm), 
it is thought that the peak position of the 11.3 fim feature may reflect the aliphatic/aromatic 
ratio. 

The 11.3 /im feature is not an isolated spectral feature. In the spectra of proto-planetary 
nebulae, weaker emi ssion features at 12.1, 12.4 a nd 13.3 /im are also seen to accompany 
the 11.3 /im feature flKwok. Volk &: Hrivnaklll999l) . It is known from laboratory studies of 
substituted aromatic molecules that the frequencie s of the OOP bending mode correlates 
with the number of adjacent H atom on each ring fjBellamvlll958h . This empirical rule has 
been translated to the correlation bet ween OOP vibrational frequenc ies and the number of 
exposed hydrogens in PAH molecules fiHudgins fc Allamandolalll999l) . In the simplest PAH 
molecule benzene, all corners are exposed and the molecule has 6 exposed H atoms. As the 
number of fused rings increases in two dimensions in a plane (growth of graphene flakes), 
the number of exposed H on each benzene ring varies. In general, the more compact a 
PAH molecule, there are fewer number of exposed H sites on each benzene ring of a PAH 
molecule. The hydrogens on these benzene units are thus labeled as solo, duo, trio, quartet, 
quintet and sextet corresponding to 1, 2, 3, 4 ,5 or 6 exposed H corners. It has been 
suggested that the wavelength of OOP mode increases monotonically from solo to sextet 
flHudgins &: Allamandolalll999[) . The observed emission features at 11.3, 12.1, 12.4 and 13.3 
/im may therefore correspond to the OOP modes of solo, duo, trio, and quartet systems. 
We should note that these OOP bending modes are also seen in more comp lex amorphous 
hydrocarbons with mixed aromatic/aliphatic structures flHerlin et al.l Il998[) . Absorption 
features at 11.36, 11.92, 12.27, and 13.27 /xm in the laboratory sample are identified as solo, 
duo, trio, and quartet OOP modes respectively. 

In addition to the narrow features mentioned above, astronomical spectra in this spec¬ 
tral region also show a broad emission plateau extending from 11 to 13/xm. This plateau 
feature has been suggested t o be due to OOP b ending modes of a collection of PAH molecules 
with non-isola ted H atoms fl Cohen et al.lll985l) . From laboratory studies of amorphous car¬ 
bon particles, iBorghesi et al.l (119871) suggested that the 11.3 /xm feature is due to combined 
contributions of both OOP aromatic CH bending and infr a red la ttice modes. Through a spa¬ 
tial and spectral study of the Orion bar, iBregman et al.l (1l989l) concluded that the narrow 
11.3/xm feature and the 11—13/xm plateau emission are produced by separate components: 
free-flying PAHs and amorphous carbon particles or PAH clusters. Highly-hydrogenated 
PAH chrysene has also been proposed as the possible cause of the 11-13 /xm broad emission 
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complex in F- and G-type stars fjJusttanont et al.lll996h. From a study of the broad emission 
plateau features observed in proto-planetary nebulae, iKwok et al.l fl200l[) suggest that the 8 
and 12-/xm plateau features are due to superpositions of in-plane and OOP bending modes 
of aliphatic side groups attached to aromatic rings. Quantum-chemical calculations indicate 
that the addition of aliphatic side groups indeed is able to produce plateau around 12 /rm 
f Sadiadi et ah 2015 ). 


In the next section, we will perform theoretical quantum chemistry calculations to de¬ 
termine the OOP bending mode frequencies of PAH molecules in the 11-12 fim region and 
explore the relationship between these modes and the astronomical 11.3 fim feature. 


3. Out-of-Plane C—H bending modes of PAH molecnles 


As the result of advances in quantum chemistry techniques, the vibrational spectrum of 
simple molecules such as PAH can be calculated very accurately. In order to explore the OOP 
bending modes of PAH molecules, we have performed quantum chemistry calculations on a 
large group of PAH molecules. A set of 60 neut ral PAH molecules are selected f rom the PAH 
on-line databases of iBoersma et al.l (120141) and iMalloci. Joblin fc MulasI (120071) . The names 
and chemical formulae of these molecules are list ed in Table 1 and their chemical stru tures 
shown in Figure [TJ We hrst us e the B3LYP (iBeckel Il993al: iHertwig fc Kochi 119971) and 
BHandHLYP hybrid functionals (lBeckelll993b() in combination with polarization consistent 
basis set PCI (jJensenll200ll. 120021) to obtain the equilibrium geometries and the fundamental 
vibrational harmonic frequencies of the molecules. The calculations were base d on density 


funct ional theory (DFT) using the Gaussian 09, Revision C.Ol software package (IFrisch et ah 


20091 ) running on the HKU grid-point supercomputer facility. The B3LYP calculations were 
done using PQS @ running on QS128-2300C-OA16 Quantu mCubeTM machine. The double 
scaling factors scheme of (iLaurv. Carlson fc WilsonI 120121) were then applied to the DFT 
vibrational frequencies. In this scheme the vibrational frequencies > 1000 cm“^ and < 
1000 cm“^ are scaled by 0.9311 and 0.9352 for BHandHLYP hybrid functionals and 0.9654 
and 0.9808 for B3LYP, respectively. The astronomical infrared emission spectra are then 
simulated by applying the Drude model at r=500 K to these scaled vibrational lines. The 
good accuracy and reliability of such DFT/Drude modeli ng in simulation of a stronomical IR 
emission bands was demonstrated in our previous study (jSadiadi et al.ll2015[) . 


All geometries have been optimized under the default criteria of the cited ab intio quan- 


^PQS version 4.0, Parallel Quantum Solutions, 2013 Green Acres Road, Fayetteville, Arkansas 72703 
URL: http://www.pqs-chem.com Email:sales@pqs-cliem.com: Parallel Quantum Solutions. 
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turn chemistry packages. The optimized geometries are all characterized as local minima, 
established by the positive values of all frequencies and their associated eigenvalues of the 
second derivative matrix. Visualizing and manipulating the results of vibrational normal 
mode analysis were performed by utilizing the Chemcr aft program. Prude model is applied 
followed by a decomposition of normal mode analysis (jSadiadi et al.ll2015l) . 


3.1. Dependence of the OOP frequencies on the number of exposed H sites 

In order to study the changes in frequency of the OOP bending modes of molecules 
with different number of exposed H sites, we have selected 7 compact PAH molecules which 
have uniform edges, i.e. each of the molecule can be labeled as having only solo, duo, 
trio or quartet sites. Because of this uniformity, these molecules are expected to have less 
complicated vibrational modes in the range of 10.8 to 15 /rm (the accepted region of spectra 
for OOP vibrations). Figure [2] shows the wavelengths of the pure OOP bending modes of 
these 7 molecules. We can see that while the duo molecules (coronene and circumbiphenyl) 
have shorter wavelength peaks than the trio (perylene and hexabenzocoronene) and quartet 
(naphthalene and triphenylene) molecules, the latter two groups are less well separated. 
While both naphthalene and triphenylene are classihed as quartet sites, they show the OOP 
bands at 12.73 and 13.44 pm respectively. Perylene, which is classihed as trio sites, shows 
the band at 12.83 pm close to that of naphthalene. 

The reason for this inconsistency can be explored by looking at the triphenylene (D3d) 
and perylene (D2h) optimized geometry and molecular symmetry. The Hs on each rings of 
these molecules are not equivalent as it is implied under the solo, duo, trio, quartet classi- 
hcation scheme. The Hs on triphenylene rings are grouped into two symmetry equivalent 
groups, i.e., 1 pair of duo H and two para solo Hs. In the same line of reasoning, the Hs 
on each of the rings on perylene are grouped into three different H atoms which are not 
symmetry equivalent of each other, although they are classihed as trio. 

For non-compact PAH molecules, there will be mixed solo, duo, trio and quartet sites. 
Our quantitative normal mode displacement vector analysis yields the number of H atoms 
with major contributions in specihc OOP vibrational mode. From these group of H atoms 
the number of solo, duo, trio and quartet H atoms are counted and converted into percentage 
values for the contribution of each edge-class in that OOP mode. The percentage of solo, 
duo, trio and quartet sites of each of the 60 PAH molecules in our sample (including both 
uniform and non-uniform edge species) is given in Table 1. 

Out of the 227 vibrational bands in the 10.8-15 pm region, there are 112 bands that are 
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found to be pure OOP vibrational modes, the other 115 are C-H OOP bands are coupled 
ring-OOP vibrations. Excluding the coupled OOP bands our vibrational analysis shows that 
all classes of H atoms could have signihcant contribution to the pure OOP bands, covering 
the entire wavelength range of 10.8 to 15 /rm. Figure [3] plots the wavelength positions of 
each of the 112 pure OOP bands, separated into solo, duo, trio, and quartet classes (shown 
in four separate panels). From the heights and concentrations of the lines, we can see that 
the solo bands are mostly clustered around 11.0 /xm, although they can extend as far as 14 
/xm. The duo modes are more evenly spread, with concentrations in particular around 11.8 
/xm, although there are also bands around 12.5 and 13.4 /xm. The trio bands are mostly 
concentrated around 13.2 /xm. The quattro modes are mostly concentrated around 13.5 /xm 
but can be seen as short as 11 /xm. 


These results does not rule out the apparent connection between the OOPs frequencies 
and the number and kind of peripheral C-H bonds. They show that correlation is not simple 
among different kinds of PAH molecule with different exposed edge. It is also even not simple 
if the number of exposed edge of a one PAH decrease by substituent. 


In our previous study (jSadiadi et al.ll2015fl on numbers of substituted ovalene structures, 
we have shown that the two OOP bands of core PAH at 11.11 and 11.78 /xm do show the 
shifts when the number of edge C-H group decreases. However the bands are shifted in 
different ways. The 11.11 /xm blue shifted to 10.98 /xm not just because of the change in the 
number of exposed edges on the ring but because of the coupling to vibartional modes of 
aliphatic groups. At the same time, the other strong OOP band at 11.78 /xm is redshifted 
instead of blue shifted when the number of exposed edge is decreased. 


3.2. PAH molecules as the carrier of the 11.3 /xm feature 

Since the OOP vibrational frequencies of PAH molecules can be calculated precisely, 
how well do they £t the astronomical 11.3 /xm feature? Our sample in Table 1 contains 60 
neutral PAH molecules containing 6 to 80 C atoms of different honeycomb structures. We 
have calculated the OOP bending mode frequencies of the PAH molecules. The simulated 
emission spectra of these molecules in the wavelength range of 11 to 11.6 /xm are plotted 
in Figure m We note that no band peaks are observed at the position of 11.30 /xm among 
this set of PAH molecules. Only two molecules in the sample, bisanthene (C 28 H 14 ) and 
Dibenzo[bc,ef]coronene (C 30 H 14 ), exhibit a band close to 11.3 /xm, at 11.38 and 11.26 /xm, 
respectively. 

Although 60 is not a large number, our results suggest that the OOP bending modes of 





neutral PAH molecules, even including PAH molecules of different geometry with solo, duo, 
trio, etc. sites, are unlikely to be able explain the observed astronomical 11.3 ^um feature. In 
the literature, support of the PAH hypothesis has cited the success in htting the astronomical 
spectra using a mixture of PAH molecules with different sizes, geometry and charged states. 
In the next section, we will explore in detail how well such httings perform in htting actual 
astronomical spectra. 


4. Fitting of the 11.3 /im feature by the PAH model 


We performed spectral htting of the observed 11.3 pm feature in hve astro nomical objects 


using the PAHdb and the IDL package AmesPAHdblDLSuite developed by iBoersma et ah 


(I2ni4h FI. Two of the objects are planetary nebulae, one is a rehection nebula, one is a post 
AGB star, and one is an active galaxy. The astronomical spectra were downloaded from 
the Spitzer Heritage Archive and the Infrared Space Observatory (ISO) Archive. We have 
subtracted the continuum using a linear ht to the line-free regions. Detail s of the htting 


proced ure can be found in iBoersma et ahl (120131. see their Sec.3). Note that iBoersma et ah 


(120131) include only pure and nitrogen-containing PAHs. In this study, all theoretical spectra 
contained in Version 2.00 of the PAHdb, including pure PAHs and heteroatom-molecules, 
was taken into account in the htting. For the spectral httings, we have employed the thermal 
approximation for all the species, in which the theoretical intensities were scaled by a Planck 
function at 500 K. A Lorentzian prohle with a band width of 15 cm“^ was assumed. 


The htting results are shown in Figures [5H9l The 5 molecular species that contribute the 
most to the htting of the feature are displayed in each of the hgures. In all these examples, 
0-containing species (in particular dihydroynaphthalene C 10 H 8 O 2 ) is a major contributor to 
the 11.3 pm feature. In some cases, its contribution is even larger than the total contribution 
from pure PAHs. Given the fact that the PAHdb contains hundreds of PAH molecules, it is 
surprising that 0-containing (and sometimes Mg- and Fe-containing) molecules are needed 
to ht the 11.3 pm feature. 


From our httings, 0-containing molecules contribute 42%, 44%, 14%, 54%, and 45% 
of the total hux of the 11.3 pm feature in BD-|-30°3639, NGG7027, NGG7023, M82, and 
HD 44179, respectively. In principle, the htting results can be used to estimate the abundance 
ratio of carbon locked in the 0-containing molecules and pure PAHs if t heir U V absorption 
cross-sections are known. The calculation method is d etailed in Tielen^ ( 2005 ) and has been 
applied to derive the abundances of Geo and PAHs (IBerne &: TielensI 120121) . However, to 


^http: //www.astrochem.org/pahdb 
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our knowledge, there is no reported measurement on the UV absorption cross-section of 
C 10 H 8 O 2 . Adopting the same value as for PAH molecules (7 x 10“^^ cm^ per C atom), we 
assume that the fraction of carbon locked in C 10 H 8 O 2 and PAHs (/) is approximately equal 
to the ratio of their contributions to the flux of the 11.3 pm feature. Using the median 
value of the flux contribution from 0-containing species (44%), we have a rough estimate of 
/ = 0.8, suggesting an 0/C ratio of about 0.1. 

We note that the Mg^"*" and Fe^^ containing molecules are similar complexes where the 
PAH is acting as vr electron donor to the metal atoms. Also, the molecules are selected by 
PAHdb to achieve a good £t and they are not necessarily physically viable molecules. For 
example, C54 (used to fit the Red Rectangle) and Cgg (used to fit M82) have very high number 
of dangling bonds and should be extremely reactive. It is unrealistic to expect molecules of 
very short lifetimes to play a role in the formation of the astronomical 11.3 pm feature in 
the interstellar medium. 


Since the code automatically selects spectra from the data pool to get the best htting 
and the data pool contains many hundreds of pure PAHs and very few heteroatom molecules, 
we can safely conclude that the inclusion of the 0-containing speci es is needed to achiev e 
a good fit. We can compar e these fits with the h ts to NGC 7023 by iBoersma et al.l (120131) . 
Examination of Figure 4 of iBoersma et al.l (120131) shows that their model spectrum contains 
features on the shoulders of the 11.3 pm feature, which are not seen in observations. This 
conhrms that the PAHdb model, even with the large number of parameters, has difficulty 
fitting the feature with pure PAH molecules and ions. 


Discussions 


A successful model for the carrier of the UIE phenomenon has to produce strong consis¬ 
tent features at 3.3, 6.2, 7.7, 8.6, and 11.3 pm, accompanied by broad plateau features around 
8 and 12 pm, and have no strong narrow features in between the strong UIE bands. The 
difficu lty in PAH molecules meeting these criteria has been discussed by ICook fc Savkallv 
(119981) . The way to get around these problems in the past is to invoke a large mix of PAH 
molecules of different masses, sizes, geometry and ionized states and £t the astronomical 
spectra by adjusting the relative abundance of hundreds of n iolecular species. Anharmoni c 
shifts and hot bands are also included to improve the hts (jPech. Joblin fc Boissell l2002l) . 
By introducing an emission model that includes a mass distribution of PAH molecules, 
temperature-depe ndent linewidths, and t he spectral energy distribution of the background 
radiation source, ICandian fc barrel (120151) are able to fit the observed prohles of the 11.3 
pm feature. For example, the asymmetric prohle of the feature is reproduced by having 
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different mass components. Recent work has expanded the PAH database t o incorporate 


dehy drogenated PAH molecnles, inclnding those with non-planar structnres flMackie et ah 


20151) . Another possibility is to include protonated PAHs. However, recent experimental IR 
spectra suggest that neither neutral-protonated (H.PAH) or cationic-protonated (H+PAH) 
forms could success fully produce the 11.3 //m feature in moderate-size PAHs such as coronene 
fiBahou et al.ll2014i) . Since astronomical UIE spectra are observed in very different radiation 
environments, it would not be reasonable to assume that the PAH population will consist of 
PAH molecules of the same mix of physical (ionization, dehydrogenation, protonation) states 
and mass distributions eyerywhere. The large number of free parameters used in these hts 
also raise doubts abou t their validity as practically any spectra can be fitted by such methods 
f Zhang &: Kwok 2015 1. 


If PAH molecules are not the carriers of the 11.3 /rm feature, what are the other possible 
candidates? Amorphous carbona ceous solids are also known to poss es s infrared vibrationa l 
based similar to the UIE bands (jpischler. Bubenzer fc Koidl 1983a b; Guillois et ah 1996 1. 
Figure [To] shows a comparison between the astronomical UIE spectrum with the laborator y 
spectrum of nanoparticles produced by laser pyrolysis of hydrocarbons flHerlin et al.lll998l) . 
The laboratory spectrum of these particles show general resemblance to the astronomical 
spectrum, not only in the major UIE bands, but also in the plateau features. The features 
are naturally broad and do not need to be artificially broadened as in gas-phase molecules. 
In addition to the 11.3 /rm feature, strong peaks are seen at 11.9 and 13.3 //m. 

In order to examine the infrared activity of hydrocarbons with mixed aromatic/aliphatic 
structures, we have calculated the simulated infrared spectra (Fig.[TT]) for 13 different exam¬ 
ples of such structures (Fig. [12]). For simplicity, we have only explored molecules that consist 
of single aromatic rings and aliphatic chains consist ing of methyl and met hylene groups. As 
such, they represent a simple subset of MAONs (iKwok fc Zhand l201l[) . The number of 
free hydrogens on all benzene rings (the numbers inside the bracket in Fig. [TTD is kept the 
same in each structure. We also include an example of pure aliphatic molecule (C 50 H 102 ) for 
comparison. 


The hrst notable difference between these molecules and PAH molecules is in the way 
that the infrared bands appear at different parts of spectra as the size of molecules increases. 
The infrared bands of PAH molecules with different sizes can appear in every part of the 
spectra as each vibrational mode has significant variations in peak wavelengths. The MAON- 
like spectra, however, have vibrational modes occurring at consistent wavelengths and are not 
highly dependent on the size of the molecules iFig. lTT]) . The 6 main features around 3.4, 6 . 2 - 
6.4, 6 .9-7.4, 9-12, 12-15, 15-20 ^m can roughly be identihed as aliphatic C-H stretching, 
aromatic C=C stretching, scissors mode of methylene C-H bonds coupled with umbrella 
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deformation mode of methyl C-H bonds, methyl and methylene deformation modes coupled 
with aromatic C-H in plane bending modes, aromatic C-H OOPs, aliphatic fragment group 
frequencies coupled with aromatic ring OOP deformation modes. However as discussed by 


Sadiadi et al.l (120151) . all of these modes have varying degree of coupling between aromatic 


and aliphatic components. The 6 peaks in Figure [TT] show some qualitative resemblance to 
the astronomical UIE spectra (lower panel. Figure [T0|) . although not exact correspondence 
in peak wavelengths. For example, the OOP bending modes of these 14 molecules range 
from 12 to 15 /rm approximately resemble the 12 /rm emission plateau feature but do not 
show a dehnite peak at 11.3 pm. Further theoretical investigation on the vibrational modes 
of MAON structures with larger sizes and complexity is needed. 


6. Conclusions 

The UIE bands are commonly found throughout the Universe and the energy emitted 
from these bands represent a signihcant output of active galaxies. Since the strengths of 
UIE bands and their distributions are used as a diagnostic tool of the physical and radiative 
environment of galactic and extragalactic objects, a correct identihcation of its origin is of 
utmost importance. Whether it originates from gas-phase molecules or complex organic 
solids can also alter our view of the chemical abundance of interstellar matter. 

The astronomical UIE bands have consistent patterns in different astronomical environ¬ 
ments and such consistency is difficult to explain by simple, gas-phase PAH molecules. The 
11.3 pm UIE feature cannot be htted by superpositions of hundreds of pure PAH molecules, 
even including PAH molecules of different sizes and charged states. In order to correctly 
identify the origin of the UIE bands, we need to explore more complex structures, specih- 
cally those with mixed aromatic/aliphatic structures. Preliminary results in this paper show 
that the MAON-like structures have consistent clusters of vibrational bands and therefore 
are viable candidates for the carrier of the astronomical UIE bands. 
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Fig. 1.— Chemical structures of the 60 PAH molecules in our sample. The bracketed 
numbers below each structure correspond to the numbers in the hrst column of Table 1. 
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Fig. 2.— Wavelengths of pure OOP bending modes of 7 simple PAH molecules with uniform 
edges. Molecules with the same number of exposed edges are plotted in the same color (duo 
in red, trio in green, quartet in blue). The benzene molecule (pure sextet, in black) is also 
plotted for comparison. The peak wavelengths of each band are also labeled on the curves. 
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Fig. 3.— The band positions of the 227 pure OOP modes of the 60 PAH molecules separated 
into peripheral H classes (from top to bottom: solo, duo, trio, and quartet). The height of 
each line corresponds to the percentage of each class. 
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Fig. 4.— Positions and profiles of the OOP bending modes of 60 PAH molecules. The 
vertical dash line represents the nominal peak position of the 11.3 pm feature. 
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Fig. 5.— Fittings of the 11.3 /im UIE feature in the planetary nebula BD+30°3639 using the 
PAHdb model. The black solid line is the continuum-subtracted ISO spectrum and the black 
dotted line is the model £t by the PAHdb. The blue line represents the total contributions 
from all the PAH ions and molecules, and the red line represents the contribution from 
oxygen-containing species. The chemical structures of the major contributors are given on 
the right side of the hgure (from top to bottom: C 10 H 8 O 2 , Cgg^, C 24 Hj* 2 ^, C 102 H 26 , and 
C 30 H 14 ). The C atoms are shown as black circles, H as small grey dots, and O atoms as blue 
circles. 
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Fig. 6.— Fittings of the 11.3 /im UIE feature in the planetary nebula NGC 7027 using the 
PAHdb model. The black solid line is the continuum-subtracted ISO spectrum and the black 
dotted line is the model fit by the PAHdb. The blue line represents the total contributions 
from all the PAH ions and molecules, and the red line represents the contribution from 
oxygen- and magnesium-containing species. The chemical structures of the major contribu¬ 
tors are given on the right side of the figure (from top to bottom: CioHg 02 , C 48 H 18 , CggH^g, 
C 47 H 17 , C 24 Hi 2 Mg+^). The C atoms are shown as black circles, H as small grey dots, and O 
and Mg atoms as dark and light blue circles, respectively. 
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Wavelength (nm) 

Fig. 7 .— Fittings of the 11.3 fim UIE feature in the reflection nebula NGC 7023 using 
the PAHdb model. The black solid line is the continuum-subtracted Spitzer IRS spectrum 
and the black dotted line is the model fit by the PAHdb. The blue line represents the 
total contributions from all the PAH ions and molecules, and the red line represents the 
contribution from oxygen- and magnesium-containing species. The chemical structures of 
the major contributors are given on the right side of the figure (from top to bottom: C 48 H 18 , 
C 10 H 8 O 2 , C 18 H 12 , C 96 H 24 Mg’'', C 32 Hi 4 '''^). The C atoms are shown as black circles, H as 
small grey dots, and O and Mg atoms as dark and light blue circles, respectively. 
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Fig. 8 .— Fittings of the 11.3 //m UIE feature in the active galaxy M82 using the PAHdb 
model. The black solid line is the continuum-subtracted ISO spectrum and the black dotted 
line is the model £t by the PAHdb. The blue line represents the total contributions from all 
the PAH ions and molecules, and the red line represents the contribution from oxygen- and 
magnesiaum-containing species. The chemical structures of the major contributors are given 
on the right side of the hgure (from top to bottom: C 10 H 8 O 2 , C 24 Hi 2 Fe'''^, Cgg^, C 54 H 18 , 
C 47 H 17 ). The C atoms are shown as black circles, H as small grey dots, O atoms as dark 
and Fe atoms as green circles. 
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Fig. 9.— Fittings of the 11.3 /im UIE feature in the reflection nebula HD 44179 using the 
PAHdb model. The black solid line is the continuum-subtracted ISO spectrum and the black 
dotted line is the model fit by the PAHdb. The blue line represents the total contributions 
from all the PAH ions and molecules, and the red line represents the contribution from 
oxygen- and magnesiaum-containing species. The chemical structures of the major contrib¬ 
utors are given on the right side of the figure (from top to bottom: C 10 H 8 O 2 , C 47 H 17 , C 54 H 18 , 
C 54 , C 24 Hi 2 Mg+^). The C atoms are shown as black circles, H as small grey dots, and O 
and Mg atoms as dark and light blue circles, respectively. 
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Fig. 10.— Comparison of the laboratory s pectrum of nanoparticles produced by laser py¬ 
rolysis of hydrocarbons (IHerlin et al.l Il998l) (top panel) with the astronomical spectrum of 
the planetary nebula IRAS 21282-1-5050 (bottom panel). The peak wavelengths (in units of 
yum) of the astronomical UIE bands are marked in red. The long wavelength (>4 /im) part 
of the spectrum of IRAS 212824 -5050 is from ISO a nd the short wavelength (<4 ^m) are 
from ground-based observations fiHrivnak et al.l 1200711 . 
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Fig. 11.— The infrared spectra of 13 simple MAON-like molecules. In spite of the different 
sizes of the molecules, the vibrational modes seem to cluster around consistent frequencies. 
The spectrum of the pure aliphatic molecule C 50 H 102 is also plotted for comparison. 
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Fig. 12.— Chemical structures of the 14 molecules used in Fig. [TTJ All the aromatic rings 
in each of the molecule has the same number of H not substituted by methyl and methylene 
groups (the number inside the brackets). One molecule (C50H102) has no aromatic rings. 
The labels isol, iso2, and iso3 refer to three different isomers of the molecule C69H84. 



Table 1 . Our sample of PAH molecules and the peripheral H sites in each molecule 


No. 

Molecule 

Formula 

solo 

duo 

trio 

percentages 

quartet 

quintet 

sextet 

1 

Benzene 

CgHe 

0 

0 

0 

0 

0 

100 

2 

Naphthalene 

CioHg 

0 

0 

0 

100 

0 

0 

3 

Anthracene 

C14H10 

20 

0 

0 

80 

0 

0 

4 

phenanthrene 

C14H10 

0 

20 

0 

80 

0 

0 

5 

Pyrene 

CieHio 

0 

40 

60 

0 

0 

0 

6 

Tertracene 

C18H12 

33 

0 

0 

67 

0 

0 

7 

Chrysene 

C18H12 

0 

33 

0 

67 

0 

0 

8 

Triphenylene 

C18H12 

0 

0 

0 

100 

0 

0 

9 

Benzo [a] anthracene 

C18H12 

17 

17 

0 

66 

0 

0 

10 

Perylene 

C20H12 

0 

0 

100 

0 

0 

0 

11 

Benzo [a] pyrene 

C20H12 

9 

33 

25 

33 

0 

0 

12 

Benzo [e] pyrene 

C20H12 

0 

17 

50 

33 

0 

0 

13 

Anthanthrene 

C22H12 

17 

33 

50 

0 

0 

0 

14 

Benzo [ghi] perylene 

C22H12 

0 

50 

50 

0 

0 

0 

15 

Pentacene 

C22H14 

43 

0 

0 

57 

0 

0 

16 

Coronene 

C24H12 

0 

100 

0 

0 

0 

0 

17 

D ibenzo [b, def ] chrysene 

C24H14 

14 

29 

0 

57 

0 

0 

18 

Dibenzo[cd,lm]perylene 

C26H14 

0 

57 

43 

0 

0 

0 

19 

Hexacene 

C26H16 

50 

0 

0 

50 

0 

0 

20 

Bisanthene 

C28H14 

14 

0 

86 

0 

0 

0 

21 

Benzo [a] coronene 

C28H14 

0 

71 

0 

29 

0 

0 

22 

Dibenzo[fg,st]pentacene 

C28II16 

13 

0 

37 

50 

0 

0 

23 

Dibenzo[bc,kl]coronene 

C30H14 

29 

29 

42 

0 

0 

0 

24 

Dibenzo[bc,ef] coronene 

C30H14 

14 

43 

43 

0 

0 

0 

25 

Naphtho[8,l?2abc]coronene 

C30H14 

7 

71 

22 

0 

0 

0 

26 

Terrylene 

C30H16 

0 

25 

75 

0 

0 

0 

27 

Ovalene 

C32H14 

14 

86 

0 

0 

0 

0 

28 

Tetrabenzocoronene 

C36H16 

25 

0 

75 

0 

0 

0 

29 

Benz [a] ovalene 

C36H16 

13 

62 

0 

25 

0 

0 

30 

Dibenzo [hi ,yz] heptacene 

C36H20 

30 

0 

30 

40 

0 

0 

31 

Circumbiphenyl 

C38H16 

0 

100 

0 

0 

0 

0 

32 

Naphth[8,2,l,abc]ovalene 

C38H16 

19 

62 

19 

0 

0 

0 

33 

Circumanthracene 

C40H16 

25 

75 

0 

0 

0 

0 

34 

Phenanthro[3,4,5,6 vuabcjovalene 

C40H16 

13 

87 

0 

0 

0 

0 

35 

Quaterrylene 

C40H20 

0 

40 

60 

0 

0 

0 

36 

Dibenz[jk,albl]octacene 

C40H22 

36 

0 

28 

36 

0 

0 

37 

Circumpyrene 

C42H16 

25 

75 

0 

0 

0 

0 

38 

Hexabenzocoronene 

C42H18 

0 

0 

100 

0 

0 

0 

39 

Honeycomb 15“ 

C46H18 

22 

56 

0 

22 

0 

0 

40 

Honeycombl6 

C48H18 

17 

66 

17 

0 

0 

0 

41 

Dicoronylene 

C48H20 

20 

80 

0 

0 

0 

0 

42 

Honeycombl? 

C50H18 

28 

55 

17 

0 

0 

0 

43 

Pentarylene 

C50H24 

0 

50 

50 

0 

0 

0 

44 

Honeycombl8 

C52H18 

33 

67 

0 

0 

0 

0 

45 

Circumcoronene 

C64H18 

33 

67 

0 

0 

0 

0 






Table 1 


Continued 


No. 

Molecule 

Formula 

solo 

duo 

trio 

percentages 

quartet 

quintet 

sextet 

46 

Honeycomb 19 

C56H20 

30 

50 

0 

20 

0 

0 

47 

Honeycomb20 

C58H20 

35 

50 

15 

0 

0 

0 

48 

Honeycomb21 

C62H22 

32 

36 

14 

18 

0 

0 

49 

Honeycomb22 

C64H22 

36 

36 

28 

0 

0 

0 

50 

Circumovalene 

C66H20 

40 

60 

0 

0 

0 

0 

51 

Honeycomb23 

C66H22 

36 

46 

0 

18 

0 

0 

52 

Honeycomb24 

C70H24 

25 

33 

25 

17 

0 

0 

53 

Honeycomb25 

C72H24 

21 

50 

12 

17 

0 

0 

54 

Honeycomb26 

C74H24 

21 

50 

12 

17 

0 

0 

55 

Honeycomb27 

C78H26 

19 

31 

35 

15 

0 

0 

56 

Honeycomb28 

C80H26 

19 

31 

35 

15 

0 

0 

57 

Honeycomb29 

C82H26 

19 

31 

35 

15 

0 

0 

58 

HoneycombSO 

C84H26 

15 

47 

23 

15 

0 

0 

59 

Honeycomb43 

C120H36 

67 

33 

0 

0 

0 

0 

60 

Kekulene 

C48H24 

50 

50 

0 

0 

0 

0 


^The number ”xx” in the label ” Honeycombxx” denotes the number of fused benzene rings in the molecular structure of the 
PAH molecule. 







